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Abstract 

^ . An analysis of the cross section for hadronic production of gluino-squark pairs 

close to threshold is presented. Within the framework of non-relativistic QCD a 
^^^^ significant enhancement compared to fixed order perturbation theory is observed 

^ ■ which originates from the characteristic remnants of the gluino-squark resonances 

\ below the nominal pair threshold. The analysis includes all colour configurations 

I of (S-wave gluino-squark pairs, i.e. triplet, sextet and 15 representation. Matching 

■ coefficients at leading order are separately evaluated for all colour configurations. 

The dominant QCD corrections, arising from initial- and final-state radiation are in- 
\ eluded. The non-relativistic dynamics of the gluino pair is solved by calculating the 

Green's function in Next-to-Leading Order (NLO). The results are applied to bench- 
mark scenarios, based on Snowmass Points and Slopes (SPS). As a consequence of 
the large decay rate of at least one of the constituents (squark or gluino), annihi- 
^ ■ lation decays of the bound state {gq) — )■ gg, g7, qZ or q'W^ are irrelevant. Thus 

the signatures of gluino-quark production below and above the nominal threshold 
are identical. Numerical results for the cross section at the Large Hadron Collider 
(LHC) at ^/s = 7TeV and 14TeV are presented. The enhancement of the total 
cross section through final state interaction amounts to roughly 3%. 



1 Introduction 



The search for new particles, predicted in supersymmetric models, is one of the most 
important tasks of the experiments at the LHC The detailed determination of particle 
masses and couplings will be crucial for the discrimination between various manifestations 
of supersymmetry (SUSY) and alternative models, even more so if one wants to distin- 
guish between the different variants of supersymmetric models, to identify the origin of 
breaking of supersymmetry and to measure the model parameters. One of the SUSY sig- 
nals will be events with missing energy or missing transverse momentum, resulting from 
cascade decays of squarks and gluinos into the lightest supersymmetric particle (LSP) 
which escapes detection. The existence of squarks, gluinos and the LSP is definitely a key 
prediction of supersymmetry. 

Depending on details of the models, in particular the mass of the LSP, gluinos and squarks 
with masses up to 3 TeV [1] could be detected. The decay rates of gluinos and squarks 
depend strongly on the ratio of squark and gluino masses. If the gluino is heavier {rrig > 
rriq), the two-body decay of the gluino into a squark and an antiquark (or its charge 
conjugate) is possible with the squark decaying into a quark plus a chargino or neutralino. 
In the opposite case (m^ < rrig) the two-body decay of the gluino is kinematically forbidden 
and the now dominant three-body decay into quark, antiquark and neutralino or chargino, 
mediated by the virtual squark, leads to a small decay rate for the gluino while the squark 
decays into the gluino plus a quark. 

The importance of squark and gluino searches has motivated a series of detailed studies 
of hadroproduction cross sections for squarks and gluinos. The lowest order has been 
evaluated long time ago [2]-|l]. Subsequently NLO SUSY-QCD corrections were calculated 
[SHT], more recently the effect of soft-gluon resummation [SHT^ was included. The present 
paper will be concerned with gluino-squark-pair production close to threshold, which 
exhibits a number of peculiar features. 

The sum of gluino and squark single decay width is expected to be much larger than the 
annihilation decay width of the squark-gluino system. In this case the bound state decay 
proceeds through the decay of the constituents and for decay rates of several GeV no well 
defined bound states exist. Nevertheless, final state interaction will lead to a significant 
lowering of the effective production threshold, to an enhancement of the cross section and 
a strong distortion of the differential cross section, in particular of the distribution in the 
invariant mass of the gluino-squark pair, with details depending on the masses and single 
decay widths of gluinos and squarks. 

This scenario has many similarities with hadronic top quark- [T514T8] and gluino-pair 
production [T9l|20] (with rriq close to rrig) close to threshold. (It differs, however, from 
gluino pair production in the case of large squark masses, where annihilation decays of (gg) 
bound states into two gluon jets could still dominate above constituent decays pTH27] .) In 
particular the distribution in the invariant mass of the gluino-squark pair can be calculated 
with similar methods. 
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A gluino-squark pair can be combined into bound states transforming under three irre- 
ducible representations, two of them attractive and one repulsive. As in Refs. p!6|[T9] 
the NLO result for the hard corrections is approximated by including the logarithmically 
enhanced terms from initial- and final-state radiation. The Green's function is evaluated 
in NLO approximation. 

The paper will be organized as follows: For a self-contained treatment we recall in Section[2] 
the quantum numbers of the bound states, discuss various SUSY scenarios and present 
the qualitative features of threshold production for the case of interest, i.e. for gluino- 
squark bound states with decay rates comparable to the level spacing of the would-be 
bound states. 

In Section [3?T] we will present the threshold enhancement (or suppression) for the vari- 
ous colour configurations using Green's functions in NLO approximation. These will be 
evaluated similar to those of the ti system discussed in Refs. [MllIT] or the gg system 
studied in [19l|20]. (For earlier investigations of squark and gluino bound state produc- 
tion at hadron colliders see [2T| - [27] .) The important difference is the non-degeneracy of 
the constituent masses, namely the gluino and squark mass. In Section [3l2] the effect of 
initial- and final-state radiation is investigated. Only S waves will be considered. The 
choice of the proper value of the strong coupling is discussed. 

Using this input, the hadronic production cross section can be evaluated in a straightfor- 
ward way in SectionHJ We limit the discussion to proton-proton collisions at 7 and 14 
TeV and give results for several of the SUSY scenarios discussed in Section|2] and com- 
pare the results to those obtained without final-state corrections. Section[5] contains our 
conclusions. 



2 SUSY scenarios, gluino-squark bound states and 
threshold behaviour 

Let us briefly recall the quantum numbers of gluino-squark pairs in the threshold region, 
classified according to their colour configurations [27]. While the combination of a spin 1/2 
and a spin particle is trivial, a colour-octet and a colour-triplet state can be combined 
into irreducible representations as follows (See also Appendix lAI): 

3®8 = 3©6©15, (1) 

while the same equation with a conjugate triplet on the left-hand side (representing an 
antisquark) results in 

3(g)8 = 3©6©T5. (2) 

The interaction can be either attractive or repulsive. In lowest order the coefficient of the 
QCD potential which governs the final state interaction is given by the expectation value 
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R 


Cr T"'^ ■ 


interaction 


3,3 


4 3 
3 2 


attractive 


6,6 


10 1 
3 2 


attractive 


15,15 


16 1 
3 2 


repulsive 



Table 1: Colour interaction of a SU{3) octet and triplet respectively antitriplet. 



of the product of the colour generators of the fundamental and the adjoint representation 
T^jFj^i. This product, in turn, can be expressed by the eigenvalues of the quadratic 
Casimir operator of the constituents, Cp = 4/3 and Ca = 3, and of the boundstate in 
representation R, Cr = (T" + F'')^^: 

ya,l . pa,2 ^ 1 ^(^a,l ^ ^a,2^2^ _ (ya,1^2 _ (^a,2)2j _ 1 (^^ _ _ ^ (3) 

The results are listed in Tab. [1] 

The interaction potential between gluino and squark in lowest order is then given by the 
"Coulomb" potential and depends on the distance r = |r| between the two constituents 



^So.W = r^, (4) 



with 



C^t^l = = 1^1,-1 for R = 3,6,15. (5) 

For the cases with negative (positive) coefficients, corresponding to attraction (repulsion), 
the cross section will be enhanced (suppressed). NLO corrections are discussed in the next 
chapter. 

The amplitudes for the production and decay of gluino and squark have to be disentangled 
according to the corresponding representation. It is convenient to define the projectors 

with the coefficients listed in Tab. [2] which obey 

i 

^ ai,bj"^bj,ck ~ '^ai,ck "-R1-R2 ) 

p£L = R- (7) 

The corresponding projectors for a general SU{N) and a sketch of their derivation are 
given in Appendix |Al together with the results for the coefficient of the potential, C^^'. 
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R 


am 






3,3 


1 

8 


3 
8 


3 
8 


6,6 


1 

4 


3 
4 


1 
4 


15,15 


5 
8 


3 
8 


5 
8 



Table 2: Coefficients of the projectors from Eq. ([6]) within SU{3). 




Figure 1: Feynman diagrams contributing at LO to gq — )■ gq. 



Restricting ourselves to the threshold region, only S-wave configurations will be retained. 
The cross section of states with higher angular momenta is suppressed by at least two 
powers of the relative velocity v of the two constituents. The classification described in 
Tab. [H is applicable to bound state and continuum production and will be important for 
the description of final state interaction. 

As mentioned in the Introduction the phenomenology of bound state production is gov- 
erned by the relative size of the sum of the decay rates of constituents, i.e. of gluino and 
squark, + F^, compared to the the level spacing AM between the ground state and 
the first radial excitation of the strongest bound colour triplet configuration. The choice 
of AM is motivated by the fact that the binding energy per se depends evidently on the 
choice of the mass definition (pole mass, potential subtracted mass, . . . ) while AM is 
convention independent. The full decay rate of the gluino-squark system is given as sum 
of the two single decay rates plus the annihilation decay rate. 

The annihilation of the gluino-squark system is dominated by the strong decay into a 
gluon and a quark. The corresponding Feynman diagrams are given in Fig. [1] We assume 
massless quarks and "Minimal Flavour Violation" (MFV) where squarks couple only the 
quarks of the same flavour. Let us in a ffist step evaluate the amplitude for gq production 
by a gluon and a quark, a result also required for the gluino-squark annihilation. At 
threshold (i.e. for pg = nig/ {nig + niq)P and Pq = niq/{nig + niq)P) the amplitude M. can 
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be decomposed into the contributions from three irreducible representations 

M {gq Mi) = iV2g^e^ (pg) u {pg, m-g) 
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Tfjl \ (P^a^n - PlU,,l) u {p„ 0) , (8a) 



u — m„ 



u{Pg,mg)!^{Pg) {PRUi^R - PM^l) u {pg,0) 



Q{mg + mq^)mg 

X 



m-g + 9m,-J pg,, + 3 (m^ + m,-J P'^,, - 3 (m^ + m,-J PgJ 



(8b) 



(Here we have used Pg ■ e{pg) = P ■ e{pg) = 0.) with the Mandelstam variables s, t, u and 
the elements of the orthogonal squark-mixing matrix 

U = f \ . (9) 

V Url Urr j 

The index i in Eqs. ( IHal) . ( ]8bl) labels the chiralities R and L. For the decay of the triplet 
and sextet ground state we find (in agreement0 with |28j ) 

^1,.., ^ _ al\Rf\^)\\mg + Qm^f 

T[[gqhr^^9q) - iumlm,{mg + m,) ' 

vm^r^-^m) = . (10) 

The elements of the squark-mixing matrix drop out after summation over quark helicities 
and in the limit of massless quarks. 

The radial part of the Schrodinger wave function at the origin can be derived from the 
one of the quarkonia or gluinonia by replacing the colour factor and the reduced mass 



with rrij-cd = mgrriq/ {jrig + nig). Numerically, for irtg ~ m ^ 1 TeV, this leads to 

a decay rate F = (75/32) a^a^ m which is of (9 (10 Me V) and thus dramatically smaller 
than the constituent decay rate discussed below. Gluino-squark annihilation into quarks 



^Wc thank Y. Kats and D. Kahawala for correspondence on this point. 



6 



Figure 2: Feynman diagrams contributing at LO to gq — ?■ 7g, Zq, W^q'. 



and electroweak gauge bosons (Fig. [2]) is further suppressed by the electromagnetic or 
weak couphng and thus completely negligible. 

The threshold behaviour, i.e. the dependence of the cross section on the invariant mass 
of the gluino-squark system is governed by the relative size of the constituent decay rate 
and the binding energy, or more precisely, the excitation energy of the dominant ground 
stat(i AM = 1 - eP with 



^[3] _ asY rrired 



2n^ - '^2) 

(Higher order corrections are included in our evaluation of the Green's function in Sec- 
tion O) 

In Tab. [3] the sum of the single decay rates, calculated with the program SDECAY [29] . 
is compared to the level spacing AM for the 17 benchmark points defined in Tab. 4 
of Ref. [20] (the spectrum of the SUSY masses has been derived with SuSpect [30]). 
Due to MFV and the structure of the proton only the gluino-squark systems of the first 
generation are of interest. For the level spacing AMgq the averaged squark mass of the 
first generation 

rUg = ^ ^, (13) 

has been used. Obviously the sum of the single decay rates is of order AMgq (benchmark 
points (a)-(i) and (o)-(q)) or far larger (benchmark points (j)-(n)). In the latter cases the 
width of the gluino is small but the one of the squark is huge in comparison with the mass 
splitting AMgq. 

The annihilation decay rate is always negligible and thus not shown here. Qualitatively 
the dominance of constituent decays is easily understood: For squark masses significantly 
larger than gluino masses the squark decay rate with its kinematically enhanced two-body 
mode becomes large. In the opposite case the gluino decay rate increases. Only in the 
degenerate case rrig ~ rriq both rates are relatively small, i.e. comparable to AM, and 
it is only in this case that resonant structures will arise. Nevertheless, even in the case 
where the decay rate is significantly larger than AM, final state interaction will lead to 



^Final state interaction for gq in the sextet representation is suppressed by powers of C'^' /C'^l = 1/3. 



benchmark 


nig 


nig 


+ r.a^ 


r- + r- 


^9 + ^d. 


^9 + 


AMgg 


point 


[GeV] 


[GeV] 


[GeV] 


[GeV] 


[GeV] 


[GeV] 


[GeV] 


(a) 


606.11 


555.39 


10.00 


5.69 


9.80 


4.83 


3.28 


(b) 


493.05 


453.62 


7.65 


4.07 


7.40 


3.38 


2.82 


(c) 


381.45 


349.95 


6.12 


3.31 


5.80 


2.78 


2.34 


(d) 


717.12 


655.32 


12.37 


7.35 


12.21 


6.34 


3.72 


(e) 


826.71 


753.94 


14.73 


9.02 


14.58 


7.85 


4.13 


(f) 


935.18 


851.53 


17.05 


10.68 


16.92 


9.36 


4.53 


(g) 


1042.60 


948.11 


19.33 


12.31 


19.22 


10.84 


4.91 


(h) 


1149.42 


1044.12 


21.57 


13.91 


21.47 


12.29 


5.29 


(i) 


936.42 


859.67 


16.42 


9.97 


16.28 


8.63 


4.55 


(j) 


802.21 


1554.00 


79.95 


65.25 


80.04 


62.71 


5.16 


(k) 


566.65 


1299.75 


79.29 


67.00 


79.40 


64.93 


4.13 


(1) 


319.59 


1055.75 


79.34 


69.31 


79.45 


67.68 


2.90 


(m) 


1030.98 


1812.80 


81.51 


64.35 


81.59 


61.35 


6.08 


(n) 


1255.61 


2073.80 


83.80 


64.20 


83.87 


60.71 


6.95 


(o) 


933.03 


841.02 


17.63 


11.36 


17.50 


10.06 


4.50 


(P) 


734.11 


753.47 


9.59 


3.34 


9.52 


2.17 


3.95 


(q) 


719.66 


664.46 


17.96 


12.85 


17.87 


11.82 


3.74 



Table 3: Comparison of the gluino mass, the averaged squark mass of the first generation, 
the sum of the single decay widths and the level spacing for 26 benchmark points defined 
in Tab. 4 of Ref. [20]. 
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a broad enhancement. Therefore results for two benchmark points will be presented in 
Section m one with Tg + Tg comparable to AM, the second one with Tg + Tg significantly 
larger. 

Hence formation of bound states [26] or of sharp resonances in the differential cross section 
(as discussed for the gluino-pair production in [IHIEO]) is not possible here. Nevertheless, 
for the benchmark points (a)-(i) and (o)-(q) the binding effects of the colour triplet and 
octet will be visible in the differential cross section. 

The technical aspects of the calculation for threshold production are very similar to the 
ones for top-quarks [T6l[I7] and gluinos [T91I20] . The cross section, differential in M, the 
invariant mass of the gluino-squark pair, can be decomposed into a factor representing the 
hard, short distance part and a factor determined by the imaginary part of the Green's 
function, evaluated at the origin. This partonic cross section is then convoluted with the 
luminosity function: 



dM ^ ' 



E 





'dCij' 


I dr 


> p 





(14) 



with 



^dOij^Ri^ /I #2 ,.2 ,.2 



dM 



s,M^/i^,/i^) 



4Im{G[^](0,M + |(r^ + r,-J)} 



red 



dr 



/ dxi dx2 fi\p{xi, fip) fj\p{x2, fJ^p) 6 {t - X1X2) . 
Jo Jo 



(15) 



The quantities S and s denote as usual the hadronic and the partonic center of mass 
energy squared, and r = s/S. The lower limit of the integration is given by p = M'^/S. 
The superscript of the Green's function refers to the colour representation R of the cor- 
responding S wave and fip and fj,R denote the factorization and renormalization scale, 
respectively. 

The only reaction present at leading order (LO) is gq — )■ gqi where in the Feynman 
diagrams of Fig. [1] initial and final states have to be exchanged. For the four squark types 
under consideration and including both squark and antisquark the hard parts of the cross 
sections are given as 



J- 



(0) _ ^[R]7r^"K/^i?) 



5(1 - z) 



(16) 



ij^R ■■^J 

with z = M'^/s. The non- vanishing normalization factors Af}^ are obtained from Eqs. ([HS 
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and (IHbl): 



•' gg 



•'^ gg 



gg 



•' gg 



gg 



•' gg 



{nig + Wiqf 
3mg{mg + rfig) 
Q{mg + mq) 
nig 

15{mg + niq) 



m-g 



(17) 



which do not depend on the squark mixing. 

The Green's function depends on the energy E = M 
decay rates. It is obtained from the non-relativistic Schrodinger equation 



Trig + niq-) and on the sum of the 



y[«](r) 



^ + ^(r, + r,-) G[^Mr,M + ^(r, + r 



5(3) 



with the potential being V^]^^^^{r), the Coulomb potential of Eqs. (jl]) and (E]). NLO 
corrections to both, the hard parts of the cross section as far as the Green's function are 
studied in the next section. 



3 Next to leading order corrections 

The NLO corrections of the cross section can be split up into those for the imaginary part 
of the Green's function ImjC'^l} and those for the hard coefficients J-'ij^r as defined in 
Eq. (USD. 



3.1 Green's function 



Following the idea of the Green's function method developed in Refs. [211122] we start 
with the interaction potential in position space. With the colour coefficients of Eq. it 
is given up to NLO 



V^"' (r) 



a. 



1 + 



a 



An 



2/3o (in + 7^,) + ai 



(19) 



The coefficient ai = 31/9Ca — 20/9Tpn/ is identical to the one relevant for the NLO 
corrections to the potentials of gluinonia and quarkonia [33]. The coefficient of the QCD- 
beta function is (3q = ll/3Cyi — 4/3Tpn/ with Ca = 3 and Tp = 1/2 and '-/e = 0.5772. The 
choice of rij = 5 active quark flavours is motivated by the scale of the Green's function 
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(see Eq. (!25|) ) for the gluino and squark masses under consideration. Effects of massive 
top quarks could be incorporated, liowever, are irrelevant at the present level of precision. 
For both LO and NLO we adopt a,(Mz = 91.1876 GeV) = 0.12018 as provided from 
MSTW2008NLO [3l] and employ the two-loop running as provided by RunDec |35| . 



The Green's function for the top-antitop system is known in compact analytic form [36] 
(see also [37]) and the result is easily applied to the present case 



Gt^] (0,i? + m^ + m,. + -(r5 + r,-J 



TT 



TT 



5'LO 



47r 



9nlo 



(20) 



with 



9lo 



fi'NLO 



(0) 



/3o 



+44F3(1,1,1,1;2,2,1-/€;1) 



+ ai 



(0) 



Kip 



(1) 



(21) 



and 



K = I 



2w 



2m 



L = In 



red 



4 mred W 



(22) 



The n-th derivative ip'^'^'^ = ip'^'^^l — k) oi the digamma function il){z) = 'jE + id/dz) lnr(2;) 
is evaluated at (1 — k). The numerical solution of the Schrodinger equation of Eq. (ITSll in- 
volves multiple poles in the binding energy which have to be resummed to single poles [36] . 
The definition of the Generalized Hypergeometric Function 4F3 together with comments 
about its numerical evaluation can be found in Ref. [SB]. In particular, attention has to 
be paid to the evaluation of 4F3 at the beginning of a branch cut for its last argument. 



3.2 Short distance corrections 

The complete NLO corrections to the hard part of the cross section involve virtual cor- 
rections to the process gq — gq plus real radiation (for gluino-pair production see [20]). 
The production of the final state gui in colour representation R, for example, will receive 
contributions from the processes gg — )■ Ru, ug — )■ Rg, uu — )■ Ru, uu — )■ Ru, uq — )■ Rq, 
uq — Rxl and qq — Rq with q^ u. 

For the present analysis the complete calculation has been replaced by an approximation 
employed in [T61IT9] for ti and gg production which includes the leading logarithms from 
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initial and final state radiation, affecting the subprocess gq gq which was present at LO 
already. The quality of this approximation can be judged e. g. from the corresponding 
analysis for top production ( [17], Eq. (14)) where both partial and complete results are 
available. In this case the size of the corrections depends on the subprocess. For gluon 
induced processes the corrections are large (35% - 40% of the corrected result) with the 
enhanced term amounting to about 3/4 of the corrections. For quark induced processes, 
however, the corrections amount to 25% - 30% with roughly equal share of enhanced and 
non-enhanced terms. One gets 



ij^R 



12 s 



SQCD 



In 



TT 



jn-g + mq) 



X l5(l-^) + ^i^^^7^[^l(^^ 



TT 



(23) 



with 



7e[^l(z) = {Ca + Cp)\2 



-c 



R 



ln(l - z) 

1 - z 

-5{l-z) 



-In 



[m-g + rUgY 



I- z 



In 



jrig 



mg)2 



(24) 



where Pq^'^^ = 3Ca — 2Tpnf is the one- loop coefficient of the SUSY QCD beta function. 
The normalization factors are given in Eq. f lT7|l . the constants Cr in Tab. [1] The 
conventional plus-distributioio has been employed to regularize singularities appearing 
at z = 1. The terms of the first line of Eq. flM|) stem from initial-state radiation and 
are therefore proportinal to the sum of the Casimirs of the adjoint and the fundamental 
representation. The first term in the second line originates from final-state radiation and 
depends on the representation of the gluino-squark state. The last term originates from 
the renormalization of the parton densities. We do not include the estimate of the virtual 
corrections based on the colour summed open production - as performed in [19] for the 
gluino-gluino case. 



4 Numerical results 



4.1 Benchmark Points 



From now on we limit our discussion to the benchmark points (p) and (q). First we focus 
on (p), where the gluino mass and the four squark masses of the first generation lie in 



^The plus-distribution follows the prescription dz 



f{z) ^ /^dzlHlfi^ [/(.)- /(I)] 
for n = 0, 1, . . . and any test function f{z). If the lower integration boundary is given by p > the 



ln"(l-z) 
1-z 



plus distribution can be replaced by the p-description via 



ln"(l-z) 

1-z 



In" 



n+1 



ln"(l-z) 
1-z 



where the latter is defined through dz 



ln"(l-z) 
1-z 



/(z)^/;dziiqi^[/(z)-/(i)]. 



12 



a range of less than 30 GeV. For the right handed squarks the sum of the decay widths 
is shghtly smaller than the level spacing AM^g, for the left handed ones approximately 
twice as big. Thus, the enhancement from the lowest lying resonance will be visible in 
the differential cross section. 

As an alternative we consider benchmark point (q). In this case the constituent decay rates 
are significantly larger than the level spacing. Final state interaction will, nevertheless, 
lead to a significant enhancement and distortion of the cross section in the threshold 
region. 

4.2 Green's Function 

As a characteristic scale for the strong coupling in the potential we have chosen the Bohr 
radius 

/xg^l = 2mred|C^t«]| (25) 

The normalized imaginary part of the Green's function at the origin is displayed in Fig. |31 
The plots involving the thresholds with up-type (down-type) squarks can be found at the 
top (at the bottom) while left-handed (right-handed) squarks can be found on the left 
(right) side. LO and NLO curves are displayed separately and the poor convergence of the 
binding energies, which had been observed for quarks [381 - HO] as well as for gluinos [26] . 
results in a shift of 2 GeV between the ground state of the most attractive colour-triplet 
configuration going from LO to NLO. The effect of the larger decay widths for the left- 
handed squarks manifests itself in the smearing of this resonance compared to the right- 
handed particles. The attractive potential of the colour-sextet representation is suppressed 
by a factor of three such that no resonances can be seen. The contribution of the repulsive 
15 configuration is negligible in the threshold region and the LO and NLO curve lie on 
top of each other. 

4.3 Cross section 

Combining the results of the previous section, the threshold behaviour of gluino-squark 
production in proton-proton collision can be computed. For the numerical evaluation of 
the strong coupling we used the program RunDec [35] and the starting values as{Mz) = 
0.13939 for the MSTW2008LO PDFs and a,(Mz) = 0.12018 for the MSTW2008NLO 
set [31]. For fiH (Eq- f EBjl ) the running is performed at the two- loop level. Decoupling 
of the top-quark at rrit = 172.0 GeV and decoupling of the gluino and squarks at is 
performed at one-loop precision (see for example [41]). 

For simplicity the squark mass entering the reduced mass of the cross section and the 
determination of the coupling is understood as averaged mass of the squarks of the first 
generation from Eq. ( IT3|) . The information of the four different squarks is only encoded 
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Figure 3: Imaginary part of the Green's function for scenario (p) for the up-type (at the 
top) and the down-type (at the bottom) squarks. LO and NLO curves are plotted for 
all three colour configurations, whereas both curves lie on top of each other for the 15 
representation. 



in the decay widths Tg. and in the total mass of the gluino-squark system entering the 
Green's function via the relative velocity v from Eq. fl22p and determining the position of 
the threshold. 

In Fig. m the LO result for the production cross section, differential in the invariant mass 
of the two supersymmetric particles, can be found. The contributions of the channels 
gq — )■ gq and gq — t- gq* are both included. We have employed the parton distribution 
functions (PDF) of MSTW2008LO [M]. Both factorization and the renormalization scale 
have been set to the hard scale 

/iF = fJ'R = ^-H = rrig + mg. (26) 

Up to NLO Green's function and hard kernel are separately renormalization scale invari- 
ant, such that the cross section involves the two scales and fin- 

The upper four plots of Fig. H] are arranged as in Fig. |3l They show the contributions 
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Figure 4: Prediction for the invariant mass distribution for scenario (p) at LO. The upper 
four plots represent the results for the configurations separately and are ordered as in 
Fig. 131 The figure at the bottom displays the individual contributions from m^, and di^^^ 
and their sum. 
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from the three representations and their sum. In the lower plot the contributions of the 
four gluino-squark channels are shown as well as their sum. It is clearly visible that, due 
to MFV and the quark content of the proton, the contributions from up-type squarks are 
twice as big as the ones of the curves involving down-type squarks. The resonances of the 
right-handed squarks nearly coincide so that a relatively sharp resonance will be visible 
in the invariant mass distribution. 

For degenerated squark and gluino masses, as it is roughly realized in this scenario (p), 
the normalization factors of Eq. f|T71) are comparable for all three colour representations. 
Therefore, the relative magnitudes of the cross sections are largely governed by the Green's 
function. Furthermore, the hard kernel does not vary strongly in in the threshold region, 
hence the shapes of the curves shown in Figs. [3] and S] exhibit a similar behaviour. 

Let us now consider the (approximate) NLO result. The NLO values for the Green's 
function will be used together with the approximation of Eqs. (1251) and fl2^ for the hard 
part of the invariant mass distribution. 

The result for benchmark point (p) is displayed in Fig. |5] for two regions of the invariant 
mass. The shapes of LO (Fig. H]) and NLO predictions (Fig. |5]) look quite similar, differing 
mainly in their relative normalization by about 20%. Furthermore, the difference between 
the masses of ul and ur, respectively, becomes irrelevant for high energies, hence the 
difference between the configurations with left- and righthanded squarks vanishes far 
above the threshold. The contribution of the qul/r final states is twice the contribution 
of the down-type squarks due to the structure of the proton. We also show the sum of the 
fixed order (FO) curves where the Green's function is replaced by its leading perturbative 
expansion in Og, furthermore, in the limit F^, F^^ — )■ 0. Hence one gets 

Im (Af + 1 (F, + F,)) } ^ ^ - (l + ^'"^ ^) • (27) 

with V as in Eq. ( 122|) . In the threshold region the four different values of the sum mg + mq. 
are clearly visible from the kinks of the dashed FO curve. Above about M = 1700 GeV the 
two curves of the FO ansatz and the one using the full Green's function roughly coincide 
and binding effects are no longer of importance. The integrated difference between the 
two functions 

/•Ml 

Aa = dM 

J Mo 

for Mo = 1415 GeV and Mi = 1650 GeV represents a measure for the effect of the final 
state interaction and amounts to a relative enhancement of about 2.57% compared to the 
total cross section of cxtot = 7.29 pb which has been calculated at NLO using Prospino2 [5]. 
It should be mentioned that also effects of squarks from the second and third generation 
are included in the total cross section. 

The sensitivity of our results on the factorization scale fip and the renormalization scale 
fiR is shown in Fig. [6], both for LO and NLO predictions. The scales ftp and fiR are 
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Figure 5: Prediction for differential cross section as function of the invariant mass for 
scenario (p) evaluated in NLO approximation as defined in Eqs. (!23l) and (l24l) for the 
four squark species. Also shown is the sum of the four squarks (solid curve) and the FO 
prediction (long dashed curve). 
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chosen to be equal and varied between /iH/2 and 2/1^^ (see Eq. [26l). As expected, the 
scale dependence is significantly reduced for the NLO approximation, amounting to less 
than iblO%. An estimate of the theoretical uncertainty may either be based on this scale 
dependence, or, more conservatively, on the difference between LO and NLO results, 
amounting to roughly 20%. 

As a second option we study the prediction for benchmark point (q) with constituent 
decay rates significantly larger than the level spacing. The Green's functions are shown 
in Fig. [7l The strong enhancement of the triplet Green's function is evident, as well as 
the strong suppression for the 15 representation. The predictions for the differential cross 
sections are shown in Fig. [8], again separately for the different squark species and the 
different representation. The step in da/dM slightly below the nominal threshold is still 
visible for the individual squark species, it is, however, smeared out when one considers 
the sum of all four channels. Nevertheless, again one observes a significant enhancement 
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of the differential cross section in the threshold region, which extends from approximately 
20 GeV below to 50 GeV above the nominal threshold. 

The NLO result for benchmark point (q) is displayed in Fig. [9] again for two regions of 
the invariant mass. The effect of the NLO correction and the dependence on factorization 
and renormalization scale is similar to the one discussed for benchmark point (p). 

5 Conclusions 

The threshold behaviour of gluino-squark production at hadron colliders has been inves- 
tigated. The NLO Green's function has been evaluated for gluino and squark in the three 
possible irreducible representations with dimensions 3, 6 and 15. The logarithmically 
enhanced part of the NLO corrections to the "hard kernel" for the evaluation of the cross 
section is included in the analysis. 

In contrast to the case of gluino pairs the constituent decay rate is generically comparable 
or larger than the binding energy, independent of the choice of the the SUSY parameters. 
Hence it is only for the exceptional case of gluino and squark masses being comparable, 
where both gluino and squark decay rates are suppressed, that pronounced resonance 
behaviour might occur. Nevertheless, attractive final state interaction, in particular for 
squark and gluino in the triplet representation, will lead to a modification of the cross 
section in the threshold region. The effective threshold is lowered by about 10 GeV from 
final state interaction and additional smearing due to the constituent decay rate. Most 
of the enhancement is concentrated in an interval of about 50 GeV around the nominal 
location of the threshold. Final state interaction, thus, leads to a significant distortion 
of the differential distribution da/dM. Compared to the total cross section evaluated 
in NLO without inclusion of the rescattering effects final state interaction increases the 
result by about 2-3%. 

A Tensor product in a general SU(N) 

The tensor product of the adjoint and the fundamental representations of SU{N) can be 
obtained following standard methods and is given by 

X N(N + 2)(N -1) N(N + 1)(N -2) , , 

{N^-1)^N = — -®— — ^^y^ -®N, (29) 

where the second representation on the right hand side of Eq. (l29l) is understood to be a 
conjugate one. (See also [27].) 

The projectors for these three representations are given as (for SU{3) see also Eq. ([6])) 
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Table 4: Coefficients of the projectors from Eq. fl30|) for a general SU{N). 



The constants a^^l, 6^^' and cl^' can be determined using the completeness relation 

E^Si = (31) 

i 

the projector properties 

'^ai,bj'^bj,ck - ^ai,ck ^RlR2 ) 

as well as the proper dimensionality of the representations 

5a.5.,PgL = dn. (33) 

The result is given in Tab. H] and the values for = 3 coincide with those of Tab. |2l 

The coefficients governing the strength of the potential, C^^\ are easily obtained by pro- 
jecting the amplitude for the one gluon exchange T^^F^^ onto the three different represen- 
tations 

dnC^^^ = PS,,,7;^K,Pii, (34) 
and the result is listed in Tab. HI 
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